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ABSTRACT 

NGC 5053 provides a rich environment to test our understanding of the complex evolution of globular 
clusters (GCs). Recent studies have found that this cluster has interesting morphological features 
beyond the typical spherical distribution of GCs, suggesting that external tidal effects have played 
an important role in its evolution and current properties. Additionally, simulations have shown 
that NGC 5053 could be a likely candidate to belong to the Sagittarius dwarf galaxy (Sgr dSph) 
stream. Using the Wisconsin-Indiana-Yale-NOAO-Hydra multi-object spectrograph, we have collected 
high quality (signal-to-noise ratio ^ 75-90), medium-resolution spectra for red giant branch stars in 
NGC 5053. Using these spectra we have measured the Fe, Ca, Ti, Ni, Ba, Na, and O abundances in 
the cluster. We measure an average cluster [Fe/H] abundance of -2.45 with a standard deviation of 
0.04 dex, making NGC 5053 one of the most metal-poor GCs in the Milky Way (MW). The [Ca/Fe], 
[Ti/Fe], and [Ba/Fe] we measure are consistent with the abundances of MW halo stars at a similar 
metallicity, with alpha-enhanced ratios and slightly depleted [Ba/Fe]. The Na and O abundances 
show the Na-0 anti-correlation found in most GCs. From our abundance analysis it appears that 
NGC 5053 is at least chemically similar to other GCs found in the MW. This does not, however, rule 
out NGC 5053 being associated with the Sgr dSph stream. 

Keywords: Galaxy: abundances - globular clusters: individual (NGG 5053) 


1. INTRODUCTION 

Globular clusters (GGs) have long served as our lab¬ 
oratories to study the dynamics and evolution of sim¬ 
ple stellar populations. With the greater availability of 
high quality and homogeneous photometric and spectro¬ 
scopic data, it has become evident that GGs are in¬ 
creasingly complex stellar populations. The Galactic 
GC population has shown that most GGs exhibit abun¬ 
dance patterns and color-magnitude diagram (GMD) 
morphology indicative of multiple popu lations (see e.g., 
iGratton et al.l I20I2I: iPiotto et al.l l2014f) . A number of 
possible sources for the gas out of which second genera¬ 
tion (SG) stars might have formed have been suggested. 
These sources include rapidly rotating massive stars, 
massive b inary stars, and intermediate-mass AGB stars 
(see e .g.. iVentura et al.l 120011: iP rantzos fc Gharbonnel 
200l iD’Ercole et al.l 120081 1201(1 120121: IdeMi^ket ah 
2009f) and all the aspects of the formation, chemical, and 
dynamical evolution are currently a matter of intensive 
investigation. 

As GGs orbit the Galaxy they are stripped of their 
stars and help populate the halo. In some instances, 
wide-held photometric surveys, such as the Sloan Digital 
Sky Survey (SDSS), have revealed large tidal tails ex¬ 
tending from GGs, giving us a snapshot of the tidal strip¬ 
ping experienced by GGs. The most famous of these ex¬ 
amples is Palomar 5 (jOdenkirchen et mi2nnih . Further 
exploration of SDSS data has revealed that other clusters 
also show large tidal features. NGC 5053 is one of these 

^ Visiting Astronomer, Kitt Peak National Observatory, Na- 
tional Optical Astronomy Observatory, which is operated by 
the Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science 
Foundation. 


clusters Isee lLauchner et aI1l2006l:Dordi fc Grebell[20T^ . 

NGC 5053 is a metal-poor GC located near the north 
Galactic cap (Z = 336°, 6 = 79°) with a Galactocentric 
radius of Roc = 17.8 kpc. Near NGC 5053, ^ 500 
pc away based on the X,Y and Z positions from Har¬ 
ris (1996, 2010 version) is NGC 5024 (M53), located at 
{I = 333°, 5 = 80°) with a Rgc = 18.4 kpc. A com¬ 
parison of the physical characteristics of these clusters is 
given in Table 1. In Figure I, a history of the previously 
determined [Fe/H] abundances of NGC 5053 are plotted 
versus the year of publication. At one point NGC 5053 
was conside red the mo st metal-poor GC in the Milky 
Way (MW) (IZinnllI985ll . 

Using SDSS photometry, iLauchner et al.l (|2006ll dis¬ 
covered a 6° tidal tail associated with NGC 5053. 
A wider held anal ysis of SDSS photometry, by 
l.lordi fc Grebell ([20101 1. confirmed the findings and fur¬ 
ther mapped the extra-t i dal fe atures of the cluster. 
A study by iChun et'al] ([201011 . using Megacam on 
the Canada-France-Hawaii Telescope, detected a tidal 
bridge-like structure between NGC 5053 and its neighbor 
M53. If this bridge-like structure is genuine, it suggests 
that the evolution of NGC 5053 has been influenced not 
only by the Galaxy, but also through interac t ions w ith 
M53. It should be noted that iJordi fc Grebell ([201011 do 
not detect this bridge-like structure. The possible com¬ 
plex dynamical history of NGC 5053 does not stop there. 
Using simulations of the tidal disruption of the Sagit¬ 
tarius dwarf galaxy (Sgr dSphl. lLaw fc Maiewskil ([2010f l 
identify NGC 5053 as a likely candidate GC that belongs 
to the Sgr dSph. 
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Table 1 

Comparison of NGC 5053 and NGC 5024 (M53) Properties 


Cluster 

R.A. 

Deck 

{l,b) 

r>a 

-^o 

jDa 

-^GC 



rb 

Ur 

c 


(J2000) 

(J2000) 


(kpc) 

(kpc) 

(pc) 

(pc) 

(pc) 

log(^) 

NOG 5053 

13:16:27.09 

-1-17:42:00.9 

336°,79° 

17.4 

17.8 

9.71 

12.39 

89.13 

0.96 

NOG 5024 

13:12:55.24 

-1-18:10:05.4 

333° ,80° 

17.9 

18.4 

2.18 

5.84 

239.9 

2.04 


Note. — “Harris (1996, Version 2010i. lMcLaughlin van der Marell ||2005D 
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Figure 1. Previous [Fe/H] determinations for NGC 5053 taken 
from the literature along with our results. Abundances determined 
from photometric data are represented by open red markers, and 
those determined by spectroscopic data are represented by solid 
blue markers. 


Motivated by the apparently complex dynamics of 
NGC 5053, and its possible association with the Sgr 
dSph, we have collected high quality spectra of red giant 
branch (RGB) stars in the cluster. With these spectra we 
have been able to directly measure the Fe, Ni, Ca, Ti, Na, 
and O abundances of clusters members in order to deter¬ 
mine if NGC 5053 exhibits the typical abundance pat¬ 
terns shown in GCs, such as the Na-0 anti-correlation. 
We will compare these abundances with those seen in 
other Galactic GCs and field stars at a similar metallicity. 
We will also present how the abundances in NGC 5053 
compare to the abundances in GCs that are generally 
considered part of th e Sgr dSph (M54, Ar p 2, Terzan 7, 
and Terzan 8; see e.g iCarretta et al.]l2014[ i. 

This paper is organized in the following manner. In 
Section 2 we will present how our observations and data 
reductions were performed, followed by a description of 
how cluster members were chosen for this study. In Sec¬ 
tion 3 we describe how the abundance analysis was car¬ 
ried out and how we assessed the errors in our measure¬ 
ments. The results of our analysis are presented in Sec¬ 
tion 4. In Section 5 we compare our results for NGC 5053 
with abundances in other MW and Sgr GCs and fields 
stars. We will conclude with possible implications of 


these abundances and future studies to come. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 

The spectroscopic data for this study were collected 
using the Wisconsin-Indiana-Yale-NOAO (WIYN) 3.5 m 
telescop43 with the Hydra multi-object spectrograph and 
a 2600 X 4000 pixel CCD (STA 1) on 2014 January 18 
and 2014 February 14. The Bench Spectrograph was 
used with the “316@63.4” grating, resulting in spectra 
with a dispersion of 0.16 A pixel'^ covering a wavelength 
range of approximately 6063 — 6380 A. The typical signal- 
to-noise (S/N) ratios of the spectra were 70-90 for the 
majority of the stars in the sample. In instances where 
a star was observed in multiple Hydra configurations, 
once on each of the nights, S/N ratios are on the order 
of 100. In addition to our program objects, high S/N 
radial velocity standards were observed at the beginning 
of each night along with a hot, rapidly rotating B star to 
use as a template for removing telluric features. 

The data were reduced using the standard IRAJ13 soft¬ 
ware to perform bias subtraction, flat fielding, dark sub¬ 
traction, and dispersion correction with ThAr compari¬ 
son lamp spectra. The full integration time (4.5 hr per 
configuration) was broken up into six separate exposures 
to reduce the effects of cosmic rays in the final combined 
spectra. The processed spectra from each exposure were 
combined using the IRAF task scombine. 

2.2. Target Selection 

The stars for this study were selected based on 
their location on t he B vs {B — V) CMD from 
iSaraiedini fc Mdonel (IT^ . as shown in Figure 2. Our 
initial sample only included those stars that lie along the 
RGB in the CMD with a.V< 15.5. At the spectral reso¬ 
lution and high S/N ratio needed to accurately measure 
abundances, stars fainter than this would require unre¬ 
alistically long exposure times. The stars that met these 
criteria, and were observed with Hydra, are represented 
by open diamonds in Figure 2. A summary of the obser¬ 
vational data for our sample is given in Table 2. There 
we list the positions, photometric values, radial veloci¬ 
ties, S/N of the spectra, and cluster membership of each 
star. The positions of the stars observed in this study are 
marked on a DSS image of NGC 5053 shown in Figure 
3. 


^ The WIYN Observatory is a joint facility of the University 
of Wisconsin-Madison, Indiana University, the National Optical 
Astronomy Observatory and the University of Missouri. 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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B-V 

Figure 2. V vs (B — V) CMP for NGC 5053. Data taken from 
ISaraiedini &: Milond dl995l') . The points overplotted with open di¬ 
amonds are those stars for which we collected spectra. These are 
the stars that met the photometric criteria described in the pa¬ 
per. The point plotted as a star marks the location of a star that 
photometrically appeared to be part of the cluster, but its radial 
velocity ruled it out as a member. 
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Figure 3. A Digitized Sky Survey (DSS) image showing 
NGC 5053. The targets for this study are highlighted by small 
circles, triangles, and a square. The triangles mark stars for which 
we were able to measure Na abundances, and the circles mark 
the stars for which we were not able to do so. The location of 
the non-member that was observed is marked with a square. The 
dashed circles represent the core, half-mass, and King truncation 
radius, respectively, going from smallest to largest. The references 
for these radii are given in Table 1. 

Radial velocities were determined for the observed 
sample using the fxcor task in the NO AO package in 
IRAF. For each star we measured its radial velocity with 
four radial velocity standards observed at the beginning 
of each night. The typical standard deviation of the 
four individual measurements for a given star was 0.3 



Figure 4. The radial velocity (RV) distribution of the stars we 
observed with Hydra. The dashed line in the plot is a normal 
distribution fit to the histogram shown. The single error bar is 
centered at the cluster average RV, and the length of the error bar 
is the standard deviation of the RV sample. The average RV value 
and standard deviation for the sample are given on top of the the 
error bar. In measuring the RVs in our sample, we found one field 
star that met the photometric criteria, but was easy ruled out as 
a cluster member with an RV of -27.1 kms~^. 

km s“^ and each measurement had an average error of 2.2 
kms“^, as estimated by fxcor. We adopted the weighted 
mean of the four measurements as the final radial veloc¬ 
ity value for each star, which is listed in Table 2. The 
error associated with each velocity in Table 2 is the error 
on the weighted mean. We also denote which stars are 
cluster members and which are non-members with a Y 
and N flag, respectively. The average radial velocity of 
cluster members is 42.0 ±1.4kms“^. Our sample aver¬ 
age radial velocity is in agreement with the previously 
determined value for the cluster Vr = 42.6 ± 0.3kms“^ 
(jKimmig et al.l[^ni4[l . The radial velocity distribution of 
the stars in our sample is shown in Figure 4. The only 
non-member we observed (Star 22), as listed in Table 2, 
is not included in Figure 4 because it sits too far away 
from the cluster distribution, and would skew the x-axis 
of the plot. 

3. ANALYSIS 

3.1. Atmospheric Parameters 

Initial estimates of the atmospheric parameters 
(Tf,ff, BCv) were determined fro m the B,V,I photom¬ 
etry (ISaraiedini fc Milond 1199511 and 2M ASS J, H, K 
for ea ch star using the relations found in lAlonso et al.l 
()1 9991 1 ■ using the (B — V), (V — K), {J — iJ), and 
{J — K) colors . The colors were d ereddened with an 
E(B-V) = 0.04 (jLauchner et al.ir2006ll . which is a consis¬ 
tent value across the literat u re, fo llowing the relations 
found by iRieke fc Lebofsk\^ (|1985ll . Surface gravities 
were calculated based on the derived effective temper¬ 
atures and bolometric corrections assuming a mass of 
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Table 2 

Observational Data 


ID 

R.A. 

J2000 

Decl. 

J2000 

V 


(B-V) 

o 

1 

o 

'Vf 

(km s“^) 

Vr 

(km s“^) 

S/N 

Member 

Y/N 

2 

13:16:26.71 

-1-17:41:01.6 

14.03 

-2.23 

1.09 

2.71 

44.2 

0.8 

83 

Y 


3 

13:16:12.23 

-1-17:46:22.9 

14.08 

-2.18 

1.08 

2.73 

44.0 

1.0 

70 

Y 


4 

13:16:20.60 

-1-17:42:46.4 

14.45 

-1.81 

0.99 

2.58 

42.2 

0.8 

82 

Y 


6 

13:16:35.95 

-1-17:41:12.8 

14.56 

-1.70 

0.97 

2.55 

42.2 

0.6 

108 

Y 


9 

13:16:26.45 

-1-17:40:07.6 

14.68 

-1.58 

0.89 

2.41 

43.8 

1.1 

108 

Y 


10 

13:16:24.50 

-1-17:39:50.0 

14.82 

-1.44 

0.90 

2.44 

40.6 

1.8 

50 

Y 


11 

13:16:57.64 

-1-17:38:05.1 

14.83 

-1.43 

0.89 

2.43 

40.5 

0.7 

71 

Y 


12 

13:16:38.19 

-1-17:40:51.8 

14.91 

-1.35 

0.95 

2.51 

42.3 

1.0 

109 

Y 


14 

13:16:37.90 

-1-17:45:35.5 

15.04 

-1.22 

0.88 

2.39 

40.1 

1.2 

80 

Y 


15 

13:16:20.74 

-1-17:41:05.9 

15.10 

-1.16 

0.85 

2.36 

43.1 

0.9 

83 

Y 


16 

13:16:25.56 

-1-17:44:02.3 

15.15 

-1.11 

0.77 

2.20 

42.6 

1.4 

79 

Y 


17 

13:16:18.18 

-1-17:39:05.7 

15.19 

-1.07 

0.78 

2.20 

42.6 

1.2 

121 

Y 


18 

13:16:32.51 

-1-17:42:17.4 

15.22 

-1.04 

0.77 

2.16 

41.7 

1.1 

60 

Y 


19 

13:16:34.20 

-1-17:43:53.1 

15.24 

-1.02 

0.83 

2.27 

40.8 

0.8 

70 

Y 


20 

13:16:22.27 

-1-17:41:53.6 

15.25 

-1.01 

0.83 

2.53 

42.1 

0.8 

82 

Y 


21 

13:16:22.13 

-1-17:38:50.1 

15.26 

-1.00 

0.83 

2.31 

40.3 

1.2 

73 

Y 


22 

13:16:28.55 

-1-17:41:12.6 

15.35 

-0.91 

0.81 

1.82 

-27.1 

0.7 

80 

N 


23 

13:16:33.91 

-1-17:42:22.9 

15.37 

-0.89 

0.80 

2.26 

39.7 

1.4 

64 

Y 




Table 3 






Table 4 





Derived Atmospheric Parameters 




Equivalent Widths 



ID 


BCv 

log(g) 


Vt 

Element 

A 

log(ff/) 

EP 


EW 


K 


cm s~^ 


kms“^ 


(A) 


(eV) 


(mA) 

2 

4460 

-0.53 

0.80 


1.67 

Fe I 

6136.61 

-1.402 

2.45 


99.9 

3 

4400 

-0.55 

0.80 


1.73 

Fe I 

6136.99 

-2.950 

2.20 


32.2 

4 

4500 

-0.48 

1.10 


1.67 

Fe I 

6137.69 

-1.402 

2.59 


999.9 

6 

4580 

-0.46 

1.10 


1.85 

Fe I 

6151.62 

-3.295 

2.18 


21.6 

9 

4660 

-0.42 

1.20 


1.85 

Fe I 

6173.33 

-2.880 

2.22 


36.1 

10 

11 

4.'S7n 

0 4fi 

1 90 


9 DO 

Note. — Table 4 is published in its entirety in the electronic edition of 

12 

44Qn 

0 !S1 

1 90 


9 OR 

the Astrophysical Journal. A portion is shown here for guidance regarding 








14 

4630 

-0.44 

1.40 


1.65 







15 

4630 

-0.43 

1.40 


2.09 







16 

4760 

-0.38 

1.50 


1.80 







17 

4750 

-0.39 

1.50 


1.41 

rameters because it was determined to not be a cluster 

18 

4850 

-0.35 

1.69 


1.34 

member from its radial velocitv. 




19 

4650 

-0.43 

1.40 


1.65 







20 

4700 

-0.33 

1.50 


2.11 







21 

4680 

-0.41 

1.50 


1.19 

3.2. Equivalent Width (EW) Measurements 


22 

91 

4700 

-0 41 

1 .^0 


1 .^n 

The iron (Fe), nickel (Ni), titanium (Ti), and calcium 


O.SMq for every star in the sample and a distance mod¬ 
ulus of (to — M)v = 16.26 (iSaraiedini et al.l[2007ll . Due 
to the limited number of measurable Fe I and Fe II lines 
in our spectral region, it was not possible to adjust the 
Teff and log(g) from their photometrically determined 
values using excitation or ionization equilibrium. 

The microturbulent velocity (vt) of ea ch star was de¬ 
termin e d using the rela ti( 3nship s foun d iniCarretta et al. 
(2004), iJohnson et al.l (1200811 . and IPilachowski et al. 
19961) . The initial estimate of Vt was an average of the 


velocities resulting from the relations found in the three 
references listed above. The Vt was later adjusted to 
minimize the dependence of derived [Fe/H] abundance 
on Fe I line strength. The final adopted atmospheric 

parameters a re listed in Table 3. _ MARCS LTE model 

atmospheres (|Gustafsson et al.l I2008D were interpolated 
to the photometrically determined atmospheric parame¬ 
ters in order to be used in the abundance measurements. 
We did not derive these stellar parameters for Star 10 be¬ 
cause its spectrum was too noisy to use for abundances 
measurements. Star 22 is also listed without these pa- 


EW measurements of individual Fe I, Ni I, Ti I and Ca I 
lines, respectively. The spectra in our sa mple were com¬ 
pared to a high resolution solar spectrum ([Wallace et al.l 
120111) to confirm the location of lines and possible blends. 
EWs were measured using the splat task in IRAF on the 
continuum-normalized spectra. The final line lists were 
complied with log{gf) values and excitation potentials 
taken from the Gafo-ESO Survey (Heiter et al. in prep). 
The EW measurements and atomic parameters for all 
lines used are available electronically. The abundances 
were determined using the abfind tas k in the 2010 Au¬ 
gust version of MOOG (ISnedenlll973D. We assumed th e 
solar abundances as found bv lAnders fc Grevess^ (|1989D . 

3.3. Spectral Synthesis 

The sodium (Na), oxygen (O), and barium (Ba) abun¬ 
dances were determined using spectral synthesis to ac¬ 
count for blends, molecular bands, and in the case of 
Ba, hyperfine structure and the isotopic mix of the ele¬ 
ment. Synthetic spectra were generated using the synth 
task in MOOG, with MARGS model atmospheres that 
matched the atmospheric parameters of each star, and an 
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Wavelength A Wavelength A 


Figure 5. An example of the spectral synthesis performed to determine the Na and O abundances for Star 6. Left panel: the spectral 
region around the Na I line near 6160 A. The strong line located at approximately 6162 A is a Ca I line. Right panel: the spectral window 
near the [O I] line near 6300 A. In each panel the data are plotted as black squares. The spectrum in the right hand panel has been 
corrected for telluric features. 




Wavelength A Wavelength A 


Figure 6. As in Figure 5, but for Star 12 with similar atmospheric parameters but very different oxygen abundance. 
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[Fe/H] abundance as determined through EW measure¬ 
ments. Assuming Gaussian broadening, the FWHM of 
the synthetic spectral lines was set to match the FWHM 
of unblended lines in the observed spectra. The final 
abundances were determined by creating synthetic spec¬ 
tra with a range of abundances for the species and line in 
question, and choosing the synthesis that best matched 
the observed spectra. An error was assigned to these 
abundances by determining how much the abundance 
could be varied before the synthetic spectrum would 
greatly deviate from the observed lines. Examples of the 
spectral synthesis are shown in Figure 5 and Figure 6 . 

In our spectral region, there are two Na I lines that we 
used to determine the [Na/Fe] abundance at 6154 A and 
6160 A. For many of the stars, the line at 6154 A was 
too weak to measure and the final [Na/Fe] abundances 
were based on the single line at 6160 A. In the cases 
where both lines were present, an average of the [Na/Fe] 
resulting from each line was used as the final abundance. 
There were also cases where both of these lines were too 
noisy to determine an [Na/Fe] abundance for the star. In 
total, we were able to determine the [Na/Fe] abundance 
in 11 of the 18 stars observed. 

Oxygen abundances were determined from the [O I] 
line at 6300 A. In this region of the spectrum there is 
a significant amount of contamination from telluric fea¬ 
tures. To remove these features, we divided the observed 
spectrum by the scaled spectrum of a hot B star, which 
acted as our template for the telluric features. The syn¬ 
thetic spectra were then compared to these telluric cor¬ 
rected spectra when determining the final [0/Fe] abun¬ 
dances. As a test, we also determined the [0/Fe] abun¬ 
dances from uncorrected spectra and there were no dif¬ 
ferences in the results. This indicated to us that the [O I] 
line was not significantly affected by the neighboring tel¬ 
luric lines. There is also a strong night sky emission line 
near this [O 1] feature. The radial velocity shift of the 
NGC 5053 stars moved the sky line off of the [O I] fea¬ 
ture. In some cases, however, the strength of the sky line 
varied enough so that the sky subtraction was not good 
enough to completely remove the line. In these cases the 
residual from the sky line was large enough to affect the 
[O I] feature. It is for that reason that we were only able 
to measure [0/Fe] abundances in 8 stars. 

The [Ba/Fe] abundances were determined from the Ba 
II line at 6141 A. In order to create an accurate syn¬ 
thetic spectrum including this feature, we had to assume 
an isotopic mix of Ba. This is complicated by the fact 
that Ba is both an s- and r-process neutron capture ele¬ 
ment. At low metallicities, such as those in NGC 5053, 
the Ba production is thought to be dominated by the r- 
process because of the slower timescales over which the 
s-process would contri bute to the overall Ba production. 
Following the work of iMcWilliaml (|1998fl . we only con¬ 
sidered the following isotopes Ba,^^® Ba, with 

a ratio of 40 : 32 : 28, respectively. We also adopt the 
r-process solar composition of log(e 0 (Ba)r) = 1-46, as 
found by iMcWilliam (1199811 . It should be noted that 
there is an Fe 1 line at 6140 A which can blend with the 
Ba 11 line, resulting in inaccurate [Ba/Fe] abundance de¬ 
terminations. At the metallicity of NGC 5053, however, 
the line is very weak so the effect is minimal. 


3.4. Uncertainties 

To quantify the errors in our abundance measurements, 
we determined how much the abundances changed with 
variations in the atmospheric parameters of a given 
model. New models were generated with one of the 
atmospheric parameters varying, while the others were 
held constant. The abundance analysis was then re¬ 
peated with these new models until we had sampled a 
sufficient grid of the parameter space to assess the errors 
induced by their individual uncertainties. A summary 
of this error analysis is given in Table 5. The analy¬ 
sis was performed on a star with the following atmo¬ 
spheric parameters over the range indicated with each 
value: = 4760 ± 100 (K), log g, = 1.50 ±0.2 (cgs), 

[M/H] = -2.5 ±0.2 (dex),ut = 1.8 ±0.2kms-h 

In Table 5 we list the CTobs. and utotai for each abun¬ 
dance measured. For abundances determined using EWs, 
(Jobs, is the standard deviation of abundances determined 
from the number of lines listed in the table. In cases 
where only one line was measured this field is left blank. 
For those abundances that were determined with spec¬ 
tral synthesis (Na, O, Ba), (Jobs, is the error assigned to 
each abundance as described in the previous section. To 
calculate atotai, we added the uncertainties from the at¬ 
mospheric parameters in quadrature. From this analysis 
we see that our observed uncertainty is consistent with 
that expected from variations in atmospheric parame¬ 
ters. 

4. RESULTS 
4.1. Fe, Ca, Ti, Ni, Ba 

In Figure 7 we plot the [Fe/H] abundance as a function 
of V, Teff, log(g), and Vt for each star to see if there are 
trends in the measured [Fe/H] values with these param¬ 
eters. The error bars in the vertical direction represent 
the standard deviation of the [Fe/H] values given by the 
Fe 1 lines in a given spectrum. We do not find any trends 
in [Fe/H] with any of these parameters. 

The individual stellar abundances for Fe, Ca, Ti, Ni, 
and Ba are given in Table 6 with the number of lines 
averaged together to reach the final abundance. In the 
second to last row of Table 6 we list the the average 
cluster abundance for Fe, Ca, Ti, Ni, and Ba with the 
number of stars used in the average and the standard 
deviation of the measurements. The final row of Table 
6 lists the assumed solar abundances used to calculate 
the abundance ratios given in this table. In Figure 8 we 
present a box- and -whisker plot of the Ca, Ti, Ni and Ba 
abundances to illustrate the median cluster abundance of 
these species, as well as their individual distributions in 
the cluster. 

4.2. Na and O 

Using spectral synthesis we were able to determine the 
Na abundance in 11 of the 18 stars in our sample, and 
the O abundance in 8 of them. There were 7 stars which 
had both Na and O abundances measured. These abun¬ 
dances are summarized in Table 7. The [Na/Fe] versus 
[Fe/H] abundance for each star and the [Na/Fe] distribu¬ 
tion for the cluster are plotted in Figure 9 and Figure 10, 
respectively. In Figure 10, the solid black line is a gener¬ 
alized histogram of the [Na/Fe] values. The generalized 
histogram was created by treating each measurement as 
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Table 5 

Uncertainties on Atmospheric Parameters 


Ion 

T^ff ±100 

(K) 

log g ± 0.20 
(cgs) 

[M/H] ±0.10 
(dex) 

vt ±0.25 
(kms“^) 

^total 

(dex) 

No. 

Lines 

<^obs 

(dex) 

Fe I 

±0.13 

±0.02 

±0.02 

=F0.03 

0.14 

9 

0.07 

[OI] 

±0.10 

±0.07 

±0.02 

±0.00 

0.12 

1 

0.20 

Na I 

±0.10 

±0.07 

±0.02 

=F0.05 

0.13 

1 

0.20 

Ca I 

±0.11 

±0.03 

±0.01 

=F0.07 

0.13 

3 

0.09 

Ti I 

±0.12 

±0.01 

±0.01 

TO.01 

0.12 

1 


Ni I 

±0.13 

±0.01 

±0.01 

TO.01 

0.13 

1 


Ba II 

±0.07 

±0.05 

±0.05 

=F0.05 

0.11 

1 

0.20 




4400 4450 4500 4550 4600 4650 4700 4750 
Teff (K) 



log g (cgs) 


1.2 1.4 1.6 1.8 

Vt (kms“^) 


2.0 2.2 


Figure 7. [Fe/H] as a function of apparent magnitude V, log(g), and vt. The vertical error bars in each of the plots represent the 
standard deviation of the [Fe/H] abundance of a given star. 


a Gaussian with a standard deviation equal to the error 
on that measurement. These individual Gaussians were 
then added together to produce the distribution shown in 
Figure 10 by the solid line. The distribution represented 
by this generalized histogram is best fit by two Gaussians 
centered at [Na/Fe] = -0.03 dex and [Na/Fe] = 0.78 dex, 
with standard deviations of 0.30 and 0.28, respectively. 
From these figures we clearly see that there is a spread 
on the order of 0.8 dex in the [Na/Fe] distribution of our 
sample, which is larger than the expected errors on the 
measurements. 

5. DISCUSSION 


5.1. Comparison With MW 

With the abundances that we determined, we can place 
NGG 5053 in the context of the chemical composition of 
stars in the MW. To make this comparison, in Figure 
11 we plot the [Ga/Fe], [Ti/Fe], and [Ba/Fe] abundances 
of NG C 5053 with a MW s ample taken from Venn et al.l 
(|2004f) . In these plots the iVeim et al.l (12004 ) sample is 
plotted as black crosses. For each species, the individ¬ 
ual cluster members of NGC 5053 are plotted as open 
squares and the cluster average for that species is plot¬ 
ted as a filled square. From these plots we can see that 
NGC 5053 lies within the envelope defined by MW stars 


































































Table 6 

Fe, Ca, Ti, Ni, and Ba Abundances in NGC 5053 


ID 


[Fe/H] (TFe n [Ca/Fe] crca 


[Ti/Fe] (TTi n [Ni/Fe] ctni n [Ba/Fe] trea 


2 

3 

4 
6 

9 

10 
11 
12 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 

Cluster Avg. 

Iog(e0) 


-2.43 

-2.43 

-2.41 

-2.33 

-2.49 

-2.44 

-2.51 

-2.47 

-2.48 

-2.50 

-2.47 

-2.45 

-2.45 

-2.43 

-2.48 

-2.45 

7.52 


0.08 

0.08 

0.07 

0.04 

0.06 

0.07 

0.08 

0.07 

0.09 

0.07 

0.07 

0.09 

0.11 

0.11 

0.07 

0.04 


17 

11 

11 

10 

12 


13 

9 

13 

9 

9 

11 

13 

13 


15 


0.60 

0.46 

0.43 

0.49 

0.44 

0.32 

0.41 

0.51 

0.32 

0.51 

0.42 

0.62 

0.43 

0.5 

0.39 

0.45 

6.36 


0.06 

0.01 

0.00 

0.04 

0.07 

0.07 

0.07 

0.14 

0.05 

0.09 

0.09 

0.06 

0.12 

0.16 

0.18 

0.08 


5 
2 
2 
4 

6 

3 

4 

5 
2 
3 

6 

2 

3 
5 

4 
15 


0.49 

0.32 

0.09 

0.26 


0.28 

0.19 

0.39 

0.32 

0.56 


0.32 

0.27 

0.43 

0.32 

4.99 


0.11 

0.14 

0.02 


0.22 

0.07 

0.20 


0.15 


0.21 

0.11 


2 

12 


-0.06 

0.07 

-0.26 

-0.07 

0.02 

0.02 

0.15 

-0.15 

0.3 


0.28 

0.15 

0.13 

0.00 

0.03 

6.25 


0.15 


0.31 

0.02 


0.01 

0.14 


1 

1 

1 

2 

13 


-0.58 

-0.33 

-0.58 

-0.58 

-0.58 

-0.58 

-0.53 

-0.78 

-0.58 

-0.43 

-0.58 

-0.43 

-0.53 

-0.53 

-0.43 

-0.54 

2.13 


0.1 

0.1 

0.2 

0.2 

0.2 

0.1 

0.2 

0.2 

0.2 

0.2 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

15 



ID 


2 

3 

4 
6 

9 

10 
11 
12 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 


Table 7 

Na and O Abundances in NGC 5053 


[Na/Fe] 


0.7 

< 0.0 

0.6 
< 0.0 

1.0 

0.6 

0.8 

0.9 

0.0 


0.9 

0.8 


Error 


0.2 


0.20 


0.2 

0.2 

0.3 

0.3 

0.2 


0.2 

0.2 


[O/Fe] 


0.2 

0.3 

0.6 

- 0.2 

0.4 


0.75 


< 0.2 

0.4 


Error 


0.2 

0.2 

0.2 

0.2 

0.2 


0.1 


0.2 


sample, we inc l ude tw o ot her measurem e nts of Ml5 from 
iPreston et al.l (1200611 and iSobeck et al.l (1201 111, and the 
abund ances for NGC 2419 as found by I Cohen &: Kirbvl 
(j2012tl . The abundance ratios from other studies have 


Figure 8. Box- and- whisker plot of the Ca, Ti, Ni and Ba abun¬ 
dances in NGC 5053. The red line is the median of each abundance, 
the bottom of each box is the first quartile, the top of each box is 
the third quartile, and the caps represent the minimum and maxi¬ 
mum values. 

at a similar metallicity. Specifically, the alpha-elements 
in our sample, [Ca/Fe] and [Ti/Fe], are enhanced by typ¬ 
ical amounts while the [Ba/Fe] is depleted. The large 
spread in the [Ti/Fe] abundances comes from the fact 
that the Til lines were relatively weak and in many cases 
only 1 or 2 lines were available in our spectral window. 


been adjusted to our assumed solar values. 

From the [Ca/Fe] panel in Figure 11, it would ap¬ 
pear that NGC 5053 is more similar to the MW field 
population than the GC sequence populated by the 
iCarretta et al.l (|2010all sample. It is difficult, however, 
to know if this offset is real or caused by the systematic 
differences between different studies. The abundances 
plotted for M15 illustrate this point well. Between the 
four measurements for M15 on this panel, [Fe/H] ranges 
from -2.3 to - 2.64 and [Ca / Fe] ra nges from 0.11 to 0.53. 
In the case of iSobeck et al.l (1201111 the difference in abun- 


On the [Ca/Fe] vs [Fe/H] plot we include the [Ca/Fe] 
abundances for other GCs in the MW. Th e open green di¬ 
amond s mark the abundances as found bv ICarretta et ahl 
(l2010ar). The aburi dances for M 68 , M15, and M30 from 
the ICarretta et ahl (| 2010 all sample are given different 
markers to highlight the locations of the three most 
metal-poor GCs in their sample. In addition to the their 


dances could arise from the evolutionary state of the stars 
being studied, but this could not account for the larger 
offsets shown in the plot. Potential systematic differ¬ 
ences in abundance analyses can arise from a multitude 
of effects, including different measurement and analysis 
techniques, adopted stellar models, stellar atmospheric 
parameters, and atomic parameters. The samples col- 
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1.4 


1.2 


1.0 


0.8 

^ 0.6 
pS 

0.4 




1 1 1 

I , 

' ' t 

t 


-0.2 I-^^^^^^-1 

-2.60 -2.55 -2.50 -2.45 -2.40 -2.35 -2.30 -2.25 

[Fe/H] 

Figure 9. The individual [Na/Fe] abundances for each star plot¬ 
ted vs their respective [Fe/H] abundance. This plot shows that 
there are two separate groups of [Na/Fe] abundances present in 
our sample. 



[Na/Fe] 


Figure 10. Generalized histogram of the [Na/Fe] abundances. 
The thick dashed line is the two-peaked Gaussian fit to this distri¬ 
bution, with the individual Gaussian components plotted as thin 
dashed lines. The residuals between the fit and the generalized 
histogram are shown in the upper left-hand corner. 


lected to make this plot did have a few Ca I and Fe I 
lines in common with our work and these had identi¬ 
cal log(g/) values, but most lines were unique to each 
study. With these considerations in mind, we can say 
that NGC 5053 exhibits alpha enhancement that is typ¬ 


ical of MW field stars and GCs at a similar metallicity 
and is among the most metal-poor GGs in the MW. 


5.2. Comparison with Sgr dSph GCs 

Throughout the literature, M54, Arp 2, Terzan 7, and 
Terzan 8 are considered to be GGs that originated from 
the Sgr dSp h, with M54 being a t the nucleus of the dwarf 
galaxy Isee llbata et alJ 11994 iDa Gosta fc ArmandrofU 
119951 : iLaw fc Maiewskill2010ll . In Figure 12 we plot the 
[Ga/Fe] of NGC 5053 along with these four clusters ver¬ 
sus their respective [Fe/H] abundances . Also plotted is 
the MW star sample from I Venn etlil] (I2004D . In this 
figure all of the abundances have been adjusted to our 
assumed solar values. As mentioned earlier and demon¬ 
strated in Figure 11, the apparent offset of some of the 
Sgr dSph GCs from the MW GC sequence could be due 
to the systematic differences between the studies. In gen¬ 
eral, the GCs considered to be members of the Sgr dSph 
do not show [Ca/Fe] values that are drastically differ¬ 
ent from either the MW GC sequence or the field stars 
at their respective metallicities, despite the systematic 
differences between studies. If NGC 5053 were to be 
confirmed as a member of the Sgr dSph it would be the 
most metal-poor cluster known to date. 

NGC 5053 would also be a relatively high-mass GC 
member of the Sgr dSph. Plotted in Figure 13 is the dis¬ 
tribution of GC masses in the MW. The locations of M54, 
Arp 2, Terzan 7, Terzan 8, Pal 12 and NGC 5053 are 
marked on the distribution. The MW GC distribution 
shown in thi s plot results from the King model masses 
reported by iMcLaughlin fc van der Mareli (1200511 . The 
individual masses marked on the distribution were cal¬ 
culated based on their respective absolute magnitudes 
reported in Harris (1996, Version 2010) and a mass- 
to-light ratio of M/L = 2. We used these calculated 
masses so that the individual clusters would be on a con- 
sistent scale. It sh ould be noted tha t rece nt work by 
iSalinas et al.l (|2012ll and iSollima et al.l (1201411 report the 
mass of Terzan 8 to be slightly higher at \ o^{M/McC\ = 
4.5 an d 4.93, respectively. The result bv ISollima et al.l 
([201411 would give Terzan 8 a mass approximately equal 
to, or slightly larger, than NGC 5053. This figure shows 
that NGC 5053 would be the second most massive GC 
associated with the Sgr dSph, after M54. 


5.3. Na-0 Anti-correlation 

As stated earlier, it was only possible to measure the 
[Na/Fe] and [0/Fe] abundances for seven of the stars 
in our sample. Despite this limited sample, NGC 5053 
still exhibits an abundance pattern in these species that 
is consistent with the Na-0 anti-correlation found in 
other clusters with larger samples, such as M54. In 
the left-hand panel of Figure 14, we plot the [Na/Fe] 
vs [0/Fe] abundances for NGC 5053 and clusters of sim- 
ilar metallicity (NGC 7099, NGC 7078, NGC 4590) from 
iCarretta et al.l ( 2009all . This panel shows that the extent 
of the Na-0 anti-correlation present in NGC 5053 is con¬ 
sistent with that found in other metal-poor MW GCs. In 
the right-hand panel of Figure 14 we plot the [Na/Fe] and 
[0/Fe] abundances for NGC 5053 and the Sgr dSph clus¬ 
ters M54, Pal 12, Terzan 7, and Terzan 8. Pal 12 is an 
additional cluster that is considered to be a likely candi- 
date associated with the Sgr dSph by iLaw fc Maiewskil 
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-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 


0.75 - + 



-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 


[Fe/H] 


[Fe/H] 


Figure 11. The Ca, Ti, and Ba abundances of the stars in our sample plotted with a Milky Way Sample taken from [Venn et all II20Q4I) . 
plotted as black crosses. The open squares in each plot represent the individual abundances for the individual stars in our sample, the solid 
re d square marks t he ave rage of these individual abundances. Left panel: the open diamonds are Milky Way GC abundances as measured 
bv ICarretta et ^ Il2010al) for Ca. Also on th i s plot, are additional measurements for M15, with references given in the legend. While M15, 
M68, and M30 are from the ICarretta et al.l (I2010al) sample, they are given different markers to highlight their values in relationship to 


NGC 5053. 



Figure 12. A comparison of the NGC 5053 [Ca/Fe] abundance 
with the GCs traditionally thought to be associated with the Sgr 
dSph: M54, Arp 2, Terzan 7, and Terzan 8. 

(j2010r) . The abundances from other studies have been 
placed on our assumed solar scale and the [Na/Fe] abun¬ 
dances have been corrected for N LTE effects following 
the models in iGratton et al.l (I1999I1 . From this panel, we 
see that Terzan 7, Terzan 8, and Pal 12 do not show 
an extended Na-0 anti-correlation and appear to be 
composed mostly, o r entirely, of a single , first genera¬ 
tion (FG) of stars. ICarretta et al.l (120141 1 suggest that 
Terzan 8 is a GC composed of mostly FG stars, with 
only 1 star in their sample appearing to belong to the 
SG. The size of our sample of NGC 505 3 stars with Na 
and O abundances is comparable to the ICarretta et al.l 
(|2014ll sample, with 7 and 6 respectively. The NGC 5053 



\og(M/M^) 

Figure 13. Distribution of GC masses in the Milky Way. The 
locations of the 4 Sgr dSph GCs are marked in the distribution with 
open points and NGC 5053 i s marked with a filled star symbol. Th e 
distribution was taken from IMcLaughlin fc van der Marell II2005I1 . 

sample, however, shows signs of a more extended Na-0 
anti-correlation and a larger fraction of stars in the SG. 
This panel also shows that the apparent FG of stars in 
NGC 5053 has [Na/Fe] abundances consistent with those 
in Pal 12, Terzan 7, and Terzan 8. 

With so few stars in the sample it is difficult to put 
constraints on how to divide the sample into the FG and 
SG stars. Based on the [Na/Fe] distribution shown in 
Figure 10, an [Na/Fe] abundance of ~ 0.4 roughly di¬ 
vides the double peaked Gaussian distribution into its 
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□ NGC 7099 A 

NGC 7078 V 

NGC 4590 ★ 

NGC 5053 

[Fe/H] -2.34 

-2.32 

-2.27 

-2.45 


5.84 

5.17 

4.91 


-1.0 -0.5 0.0 0.5 1.0 

lO/Fe] 



o M 54 i Pal 12 NGC 5053 

♦ Ter 8 ■ Ter 7 

-1.0 -0.5 0.0 0.5 1.0 

[O/Fe] 


Figure 14. Left panel: plot of the [Na/Fe] vs. [O/Fe] abundances in NGC 5053, NGC 7099, NGC 7078, NGC 4590. The [Fe/H] abundance 
and mass for each cluster is given in the plot legend. Right panel: plot of [Na/Fe] vs. [O/Fe] for NGC 505 3 and other Sg r dSph candidate 
GCs. The abundances for M54, P al 12, Terzan 7, and Terzan 8 were taken from fUarretta et al.l l|2010bl'l . ICohenI l|200Jl . fSbordone et al.l 
II2007I) . and ICarretta et ahl II2014IL respectively. All of the abundances have been put on our abundance scale and corrected for NLTE 
affects. 


two components, with an [Na/Fe] < 0.4 being the FG 
and [Na/Fe] > 0.4 being the SG of stars. Following 
this separation criterion, our sample would contain 3 FG 
stars and 8 SG stars. The fractions of FG and SG stars 
in NGG 5 053 are consistent wit h those found in other 
MW GCs (|Carretta et aLll^OObH . In Figure 3, we show 
the distribution of these stars on the cluster field; the FG 
stars are marked with blue triangles pointing up, and the 
SG stars are marked with red triangles pointing down. 
Unfortunately, there were too few stars to explore the ra¬ 
dial distribution of the different generations in the clus¬ 
ter. Out of the 5 GCs plotted in Figure 14, NGC 5053 
is the second most massive, after M54. In this less mas¬ 
sive Sgr dSph candidate GC sample, NGC 5053 is the 
only GC to show signatures of SG stars. While this does 
not put any constraints on the likelihood that NGC 5053 
is a Sgr dSph cluster, it does illustrate how cluster mass 
affects and balanc es the presence o f multip le generations. 

A study by iSmoliLki et al.l (1201111 found that 
NGC 5053 cluster members have an approximately 
unimodal CN band strength distribution, suggesting 
only FG stars are in the cluster. While their result my 
seem to contradict what we fi nd in NGC 505 3 , this 
is n ot the case. A s note d by iSmolinski et al.l (1201111 
and iShetrone et al.l (|2010ll . their analysis of CN and 
CH molecular bands is limited by the [Fe/H] of the 
clusters. At low metallicities, such as those found in 
NGC 5053, the variations in these molecular band 
strengths will become difficult to detect, causing clusters 
to appear to be chemically homogeneous in light element 
abundances. This leads them to conclude that a low 
metallicity cluster showing a unimodal band strength 
distribution could still have chemical inhomogeneities 
and multiple populations in the cluster. 


6. CONCLUSIONS 

NGC 5053 is possibly a very dynamically complex GC 
in our Galaxy. Not only has this cluster been shown to 
have significant tidal tails, but one study has also de¬ 
tected a tidal bridge with the neighboring cluster M53. 
The already rich dynamical environment of NGC 5053 
is further complicated by its possible association with 
the Sgr dSph. The goal of this study was to determine 
if NGC 5053 shows the typical chemical signatures of 
MW GCs by collecting high quality spectra for a sample 
of RGB stars in the cluster. By using both EW mea¬ 
surements and spectral synthesis, we have been able to 
determine the Fe, Ca, Ti, Ni, Ba, Na, and O abundances 
in NGC 5053. 

The abundance results are summarized as follows: 

1. The cluster average [Fe/H] = -2.45 from 15 stars, 
with a = 0.04, making NGC 5053 one of the most 
metal-poor GCs in the Galaxy. 

2. From our target selection we determine the radial 
velocity of the cluster to be Vr = 42.0 ±1.4 kms“^. 

3. At this metallicity, the [Ca/Fe], [Ti/Fe], and 
[Ba/Fe] abundances are consistent with the abun¬ 
dances found for halo stars in the MW. Specifi¬ 
cally, the Qf-elements are enhanced and the Ba is 
depleted, relative to solar. 

4. The Na and O abundances exhibit star-to-star 
variations within our sample. These variations 
are larger than the typical uncertainty in their 
respective measurements. The [Na/Fe] distribu¬ 
tion is well fit by two Gaussians with mean val¬ 
ues separated by 0.8 dex. The Na and O abun¬ 
dances exhibit a distribution that is consistent with 
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other GCs in the MW. From these abundances, 
NGC 5053 appears to have both a FG and a SG 
of stars still present in the cluster. The relative 
fraction of FG and SG stars in our sample for 
NGC 5053 is consistent with other GCs in the 
MW. The extent of the Na-0 anti-correlation in 
NGC 5053 is similar to what is seen in other metal- 
poor MW GCs, specifically M15, M68, and M30. 

5. Like the other Sgr dSph candidate GCs, NGC 5053 
appears chemically similar (e.g., alpha-enhanced) 
to the MW field population at a similar metal- 
licity. Unlike MW GCs, however, the other Sgr 
dSph GC candidates do not show an extended Na- 
O anti-correlation or clear signs of multiple pop¬ 
ulations, apart from M54. Out of the candidate 
GCs with abundances we pulled from the litera¬ 
ture, NGC 5053 is the most massive, after M54, 
and shows the signature of the most extended Na- 
O anti-correlation. In the case of the other Sgr 
dSph candidate GCs, this could be due to the small 
sample sizes available and the influence of their rel¬ 
atively low masses. 

6. Based on these pieces of evidence we are unable 
to put further constraints on the possibility that 
NGC 5053 is a Sgr dSph cluster, so it will have to 
remain a maybe. While the Na-0 anti-correlation 
we see is similar to those seen in other MW GCs, 
and not in the low mass Sgr dSph GCs candidates, 
this does not rule out its membership with the Sgr 
dSph. This could be the result of NGC 5053 be¬ 
ing massive enough to support multiple generations 
while the less massive Sgr dSph clusters are not. 
These results do suggest, however, that the possi¬ 
ble complex dynamical history of NGC 5053 has 
not affected its ability to produce the abundance 
trends we expect to see in GCs and the MW. 

This work is the first piece of a comprehensive study of 
NGC 5053 and its neighbor M 53. Using WIYN, we have 
collected deep, wide-held photometry of the area around 
the two clusters, as well as the area between them, to 
fully charcterize the morphological properties of their 
outer regions and put further constraints on their tidal 
tails and the possible tidal bridge. Additionally, we have 
collected Hydra data for M53 that will allow us to explore 
its abundances in a similar way to this study. We will 
pair these observational data with N-body simulations of 
the clusters to create a more dynamically complete pic¬ 
ture of the evolution of these clusters and the origin of 
their tidal features. 

We would like to thank to Anna Lisa Varri for the ini¬ 
tial discussions that helped seed this project, Maria J. 
Cordero who assisted with the observations and data re¬ 
duction, and Jamie Overbeek who gave guidance on how 
to determine Ba abundances. EV acknowledges support 
by grant NASA NNX13AF45G. This material is based 
on work supported by the National Science Foundation 
Graduate Research Fellowship Program under Grant No. 
DGE-1342962 to OB. Any opinion, findings, and conclu¬ 
sions or recommendations expressed in this material are 


those of the author (s) and do not necessarily reflect the 
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